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Miniature laser sources with on-demand beam features are desirable devices for a broad range of photonic 
applications. Lasing based on direct-pump of miniaturized waveguiding active structures offers a low-cost 
but intriguing solution for compact light-emitting devices. In this work, we demonstrate a novel family of 
three dimensional (3D) photonic microstructures monolithically integrated in a Nd:YAG laser crystal wafer. 
They are produced by the femtosecond laser writing, capable of simultaneous light waveguiding and beam 
manipulation. In these guiding systems, tailoring of laser modes by both passive/active beam splitting and 
ring-shaped transformation are achieved by an appropriate design of refractive index patterns. Integration 
of graphene thin-layer as saturable absorber in the 3D laser structures allows for efficient passive 
Q-switching of tailored laser radiations which may enable miniature waveguiding lasers for broader 
applications. Our results pave a way to construct complex integrated passive and active laser circuits in 
dielectric crystals by using femtosecond laser written monolithic photonic chips. 



Miniature photonic devices, such as beam splitters, directional couplers, and mirroring modulators, have 
received a broad range of applications, including but not limited to optical telecommunications, 
quantum computing, biophotonic sensing, and information processing 17 . These devices are constructed 
based on the optical waveguide technology. Versatile guiding structures with diverse geometries have become 
basic components of modern integrated photonic circuits. Compact, integratable lasers with excellent perform- 
ance are particularly desirable as light sources for intricate photonic chips so as to reduce the scale of the devices 8 . 
Waveguide lasers have been realized in large numbers of active waveguide systems, covering broad spectral region 
ranging from visible to mid-infrared, with continuous-wave (CW) and pulsed regimes (Q-switched or mode- 
locked) 913 . Some excellent performances of obtained waveguide lasers have shown the potential as light sources 
for the development of highly integrated photonic circuits. Nevertheless, the light-beam profiles of the waveguide 
lasers cannot be manipulated in a flexible manner. And in most cases, researchers manufacture waveguide lasers 
in straight channels or planar layers that strongly limit the actual control of the spatial properties of the generated 
laser radiations 8 . In practice, versatile requirements on the laser beams are desired, in addition to the easily 
achievable single Gaussian-like profiles. 

Femtosecond laser pulses with high intensities induce extremely localized modifications of material matrix 
through strong-field ionization processes. The femtosecond laser writing has become a powerful, widely applic- 
able technique to implement one-step 3D micro-engineering of materials for a great variety of applications 14-15 , 
owing to many advantages, such as wide applicability of materials, negligible thermal-diffusion effect, and ability 
for maskless 3D processing 15 . Since the pioneering work of Davis et al. in 1996, in which optical waveguides were 
fabricated in a few glasses 17 , many research groups have utilized the femtosecond-laser writing to manufacture 
guiding devices in transparent optical materials, including glasses, dielectric crystals, ceramics, and organic 
materials, and realized applications in diverse disciplines 18 22 . Unlike the modification of glasses, in which the 
femtosecond laser pulses usually induce positive changes of the refractive index (An > 0), negative index changes 
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Figure 1 | Photonic structure characterization, (a), The fabrication schematic plots of a typical photonic microstructure with a guiding core surrounded 
by a hexagonal track array, (b), Cross-sectional sketch of the photonic microstructure. (c), Microscopic image of the cross section of the fs laser written 
photonic Element-1 in Nd:YAG crystal, (d), Measured near-field modal profiles along both TE and TM polarizations at 1064 nm. (e), All-angle light 
transmission along the transverse plane at 1064 nm. (f), Micro-photoluminescence (nPL) mappings (spatial distributions of emitted intensity, spectral 
shift, and linewidth) of the photonic structure, (g), Simulated modal profiles of the guiding microstructure along both TE and TM polarizations. 



(An < 0) are typically induced in crystalline media through the 
femtosecond laser modifications. This, to some extent, blocks the 
applicability of 3D micromachining of waveguides in dielectric crys- 
tals by using femtosecond laser writing. People have produced wave- 
guides with Type I (direct- write waveguides, guiding core located in 
the femtosecond-laser-written tracks of An > 0), Type II (stress- 
induced ones, guiding core typically located in the center region of 
two parallel tracks of An < 0), and Type III (depressed cladding 
structures, core located in the region surrounded by a number of 
tracks of An < 0) configurations in crystals 19 . Particularly, femto- 
second laser writing has been employed to successfully implement 
curved waveguiding structures and waveguide lasers in Yb:YAG 
crystals recently, which paves a way to realize lasing in non-straight 
channels 23 . In this work, we present a novel family of photonic struc- 
tures consisting of hexagonal lattices of tracks (An < 0) with induced 
defects capable of efficient light confinement in a wide spectral range. 
By changing the transverse profile of the waveguide (i.e., introducing 
new axial defects in the appropriate positions of the lattice), the light 
beam manipulation of 1 X 2 beam splitting and ring-shaped trans- 
formation are achieved. We demonstrate here that the flexibility and 
convenience of femtosecond laser micro-modification enable the 
integration of diverse elements in a single crystal wafer, providing 
the possibility of simultaneous generation and manipulation of 
guided light beams. By using direct-pump of the monolithic inte- 
grated structures, waveguide lasing has been obtained in both CW 
and pulsed (passively Q-switched by graphene) regimes. 



Results 

The structures are fabricated in the well-known Nd:YAG (neo- 
dymium doped yttrium aluminum garnet) laser crystal, which has 
a cubic lattice geometry for laser oscillations at 1.06 um. The 
femtosecond-laser pulses create micro -modifications in the irradiated 
areas of the crystal, which can be controlled by the parameters such as 
pulse energy, scanning speed, beam intensity distribution, and beam 
polarization. Figure la and lb show the schematic plots of a typical 
photonic microstructure (with a guiding core surrounded by a hex- 
agonally arrayed track lattice, hereafter referred to Element 1) and its 
cross-sectional topography, respectively. The laser inscribed tracks 
have a transverse length of 10 um, and the separation between two 
adjacent tracks is set to 10 um. The central core is approximately a 
region with area of 30 X 30 um 2 to maintain a nearly symmetric 
geometry. Figure lc depicts the optical microscope image of the cross 
section of Element 1 fabricated inside a 10-mm long Nd:YAG crystal. 
The guidance of the element is experimentally investigated by end- 
face coupling arrangement at wavelength of 1064 nm with a linearly- 
polarized CW laser. As shown in Fig. Id for the measured near-field 
modal profiles, the structure supports well-confined modes along 
both TE and TM polarizations. This feature is similar to the 
depressed-cladding waveguides in Nd:YAG 24 produced by the 
femtosecond-laser writing but significantly different from the Type 
II structures which only guide TM modes 25 . As to obtain thorough 
information of the polarization effects of the guidance, we measure 
the all-angle light transmission along the transverse plane (Fig. le). It 
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Figure 2 | Photonic structure characterization for beam splitting in passive regime, (a), Schematic images of Element 2 (top) and 3 (below), 
(b), Microscopic pictures of Element 2 (top) and 3 (bottom), (c), The prototype of the beam splitter with connection of Element 1, 2 and 3, with each 
element 3.3-mm long, (d), Simulated beam profile evolution as the 1064-nm light propagates along the photonic structure, (e), Measured intensity 
distributions of the beam splitter in both passive and active regimes at 1064 nm. 



is found that the guidance is polarization-independent (i.e., the guid- 
ance exists for the input light at any transverse polarizations), reveal- 
ing a perfect isotropic light confinement capability, which is even 
better than the normal Type III Nd:YAG cladding waveguides. It 
should be also noted that the transmitted light through the guiding 
structure preserves the original polarization. Finally, the propagation 
loss of Element 1 is determined to be —0.7 dB/cm, showing an excel- 
lent light transmission property. This value is comparable with the 
Type II Nd:YAG dual-line waveguides (0.6-1.6 dB/cm) 25 , and may be 
further reduced by increasing the guiding core volume 26 . In addition, 
the amount of the track lattice layers is an important factor for light 
field restriction: multiple-layer track lattice offers better light confine- 
ment than single-layer structure (see Supplementary Fig. 1). 

Figure If shows the confocal micro-photoluminescence ((J.-PL) 
maps (spatial distributions of emitted intensity, spectral shift, and 
linewidth) of Element 1 from the cross-sectional view. These fluor- 
escence maps are obtained by analyzing the spatial variation of the 
spectral properties of the fluorescence line at around 940 nm, which 
belongs to the Nd 3+ ion transition band of 4 F 3/2 — > 4 Ig/2, which has 
been demonstrated to be a versatile fluorescence probe for the detec- 
tion of slight structural modifications in the YAG network. In the 
femtosecond-laser induced filaments, an obvious modification of PL 
emission is observed, however, inside the central region (i.e., wave- 
guiding core), the fluorescence features are well preserved with 
respect to the bulk. This confirms that the photonic structure 
Element 1 is an effective gain channel for the optical amplification 
and lasing. In addition, based on the uPL maps, the refractive index 
change in the tracks (the contrast between individual modification 
and the unprocessed crystal) is determined to be An ~ -4 X 10~ 3 , 
considering the combination of the damage induced refractive index 
reduction and the stress induced refractive index increment at the 
damage tracks and their surroundings 27 . With the reconstructed spa- 
tial distribution of the refractive index, the modal profiles (Fig. lg) of 
Element 1 along both TE and TM polarizations are calculated by 
Rsoft® software based on finite-difference beam propagation method 
(FD-BPM), showing reasonable agreement with the measured modal 
intensity distributions (Fig. Id). 



The light propagates in a straight way through Element 1. In order 
to implement beam manipulation, specially designed photonic ele- 
ments are required to connect the basic one (i.e., Element 1) in a 
reasonable manner. Since the light propagates along the central core, 
one can introduce additional axial "defect" lines in appropriate posi- 
tions to steer the beam propagation and tailor the intensity beam 
distribution. In addition, the area of the central-core region may be 
also changed to allow smooth and effective confinement of light 
beams. In this work, we present the 1X2 beam splitting and ring- 
shaped transformation, respectively, through two configurations of 
monolithic photonic structures. 

Figures 2a and 2b show the schematic (left) and microscopic 
(right) of the Elements 2 and 3, which have been fabricated by intro- 
ducing defect tracks in the core region of Element 1. The prototype 
with combination of Elements 1, 2 and 3, which are connected in 
sequence, is specially designed for 1X2 beam splitting (Fig. 2c). 
Figure 2d shows the beam profile evolution (simulated by FD-BPM) 
of 1064-nm light propagating along the three-element photonic 
structure. With Element 1, 2 and 3 connected in the sequence indi- 
cated by the scheme, a clear beam-profile splitting is observed at the 
waveguide output (Fig. 2e for experimental results) with measured 
intensity splitting ratio of 0.505 to 0.495 for the two arms, showing 
excellent performance for a Y-branch-like function. 

Similarly, we perform the beam manipulation from a single-spot 
Gaussian profile to a ring-shaped one through connection of 3 ele- 
ments. Figures 3a and 3b show the schematic (left) and optical image 
(right) corresponding to the Elements 4 and 5. The combination of 
Elements 1, 4, and 5 in a joint structure (see Fig. 3c for the prototype 
scheme) is designed for ring-shaped beam transformation. Figure 3d 
depicts the simulated beam profile evolution of light propagation 
through the structure at 1064 nm. It could be found that with the 
combination of Elements 1 and 4, a ring-shaped profile starts to be 
established; after additional connection of Element 5, the beam shape 
becomes a nearly perfect ring (Fig. 3e for experimental data). Please 
note these results are obtained experimentally in the "passive" opera- 
tion regime, i.e., under input and output signals of a 1064-nm laser 
beam. In addition, it is worth pointing out that the output powers of 
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Figure 3 | Photonic structure characterization for ring-shaped transformation in passive regime, (a), Schematic images of Element 4 (top) and 5 
(bottom), (b), Microscopic pictures of Element 4 (top) and 5 (below), (c), The prototype of the ring-shaped photonic structure with connection of 
Element 1 , 4 and 5 with each element 3.3-mm long, (d), Simulated beam profile evolution as the 1064-nm light propagates along the ring-shaped photonic 
structure, (e), Measured intensity distributions of the ring-shaped structure in both passive and active regimes at 1064 nm. 



the single-element waveguide and the three-element prototypes (i.e. 
1 + 2 + 3 or 1 + 4 + 5) are nearly the same, showing that, in passive 
regime, the axial defects do not introduce additional losses, and the 
beam manipulation can be implemented in the photonic structures 
efficiently (see Supplementary Fig. 2). 

In practice, for more compact geometry of photonic devices, the 
direct-pump lasing is on demand from waveguiding structures, i.e., 
the systems operate in the "active" regime. In straight channel wave- 
guides, efficient lasers have been realized in many gain materials 8 ; on 
the contrary, lasing structures capable of beam generation and 
manipulation were reported in limited literatures, such as curved 
waveguides in Yb:YAG crystals 23 . In this work, we apply an optical 
direct-pump to the monolithic combination of the photonic ele- 
ments to achieve manipulated guided lasers in CW regime. 
Figure 4a shows the schematic of the waveguide lasing experiment 
setup. The pump beam is taken at a wavelength of 808 nm from a 
ThSapphire CW laser. The output lasing is at the wavelength of 
1064 nm, which corresponds to the main emission line of the 
Nd 3+ ion transition band 4 F 3/2 — » 4 In/2- The obtained modal profiles 
at 1064 nm in the active regime are as same as those achieved in 
passive regime at same wavelength for corresponding combinations 
of photonic elements. The measured CW laser spectrum is depicted 
in Fig. 4b, centering at 1064 nm. Figure 4c shows the output laser 
powers at 1064 nm as a function of launched pump power for the 
three-element beam splitting and ring-shape-transformation struc- 
tures, respectively. For comparison, the lasing performance curve for 
a single Element 1 structure is also depicted in this figure. As the 
launched pump power is above a threshold (180 mW), lasing oscilla- 
tions are observed in the photonic structure systems. For the beam- 
splitted structure (Element 1 + 2 + 3), the laser shows a slope 
efficiency of 32% and a maximum output power of 241 mW. For 
the ring-shaped structure (Element 1 + 4 + 5), the laser operates 
with a slope efficiency of 29%, climbing to a maximum output power 
of 224 mW in case of a pump power of 0.95 W. Compared with the 
lasing performance based on a single Element 1 (with slope efficiency 
of 45% and maximum output power of 350 mW), i.e., a straight 
photonic structure, there are 30% and 35% loss of maximum output 



power for the beam splitted and ring-shaped configurations, respect- 
ively. This should be mainly attributed to the additional losses from 
the splitted lasing cavity, i.e., non-perfect Fabry-Perot cavity, as the 
losses from the introduction of axial defects are negligible in passive 
regime. These loss values in active regime are comparable to those of 
the curved Yb:YAG waveguide lasers (with a loss of —20% with 
respect to the straight waveguide lasers) 23 . Nevertheless, such losses 
could be further reduced by optimizing the parameters of the 
femtosecond-laser writing and dimensions of each element. In addi- 
tion, the lasing performance could be also improved by an optimized 
pump system, reaching higher slope efficiency and Watt-level max- 
imum output powers. 

Pulsed lasers are required in high-intensity applications, such as 
nonlinear optics. Graphene is a low-cost 2D nano-scale material, 
which has been widely applied in many fields 28,29 . It has been proved 
to be an ideal material for fabrication of ultra-broadband saturable 
absorbers (SAs) owing to the linear energy dispersion relation and 
fast responses 30,31 . In this work, a graphene SA mirror is attached on 
the crystal end face to obtain passive Q-switched lasers from the 
photonic structures. The measured Q-switched laser spectrum and 
the pulse train are depicted in Fig. 4b. The central wavelength of 
1064 nm is consistent with the case of CW operation. The inset in 
Fig. 4b shows the pulse train of the Q-switched waveguide laser. 
Figure 4d illustrates the average output powers of the single 
Element 1, 3-element beam splitting, and ring-shape-transformation 
structures, as functions of launched pump power at 808 nm, respect- 
ively. From the linear fit of the experimental data, we have deter- 
mined that the lasing thresholds for pulsed lasers are all around 
265 mW, which is slightly higher than those of CW regimes. 
Above the lasing threshold, the Q-switching operates pulsed lasers 
with slope efficiencies of 35%, 24%, and 22% and maximum output 
powers of 243, 163, and 148 mW for these three configurations, 
respectively. In addition, by varying the pumping power, the repe- 
tition rates of pulsed lasers from the single Element 1, 3-element 
beam splitting, and ring-shape-transformation structures could be 
tuned from less than 1 MHz to above 4 MHz (as depicted in Fig. 4e), 
when the pumping power is increased from 265 mW to 950 mW. 
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Figure 4 | Lasing with beam manipulation, (a), Schematic of experimental arrangement for the laser generation in Nd:YAG photonic structure in both 
CW and pulsed regimes, (b), Laser emission spectrum from the Nd:YAG waveguide in CW and graphene-SA Q-switched regimes. The inset shows the 
pulse train of the Q-switched waveguide laser, (c), Output laser powers from Element 1, beam splitter and ring-shaped structures as functions of the 
launched pump in CW regime, (d), Output laser powers from Element 1, beam splitter and ring-shaped structures as functions of the launched pump in 
pulsed regime, (e), Repetition rate of the Q-switched pulsed lasers, (f), Q-switched pulse durations as functions of launched pump power, (g), Q-switched 
single-pulse energies as functions of launched pump power. 



Meanwhile, the pulse duration varies from 300 ns to around 70 ns 
(Fig. 4f). Finally, the maximum single pulse energy is —40 nj for the 
3-element beam splitting and ring-shape-transformation structures, 
which is lower than that (—55 nj) for the single Element 1 (Fig. 4g). 

Discussion 

We have designed and fabricated in a Nd:YAG laser crystal wafer a 
novel family of monolithic photonic structures consisting of hexa- 
gonal lattices of laser-damage bits with induced axial defects that act 
as guiding cores. The integration of a few specially designed photonic 
elements with defect lines has been demonstrated to be an efficient 
and simple way to build compact structures capable of efficient beam 
profile manipulation, including basic functions of beam splitting and 
ring-shaped transformation. Efficient lasing in both CW and pulsed 
regimes are implemented, demonstrating that the monolithic photo- 
nic structures are capable of not only light manipulation but also 
efficient light generation. Our work paves a solid way for a wide 
variety of applications in many disciplines. On one side, from the 
fabrication point of view, our work shows the capability of femto- 
second laser microfabrication to produce more complex waveguid- 
ing devices by integration of more designable elements in a single 
crystal chip, which would further enlarge the scope of promising 
applications related to this fabrication technique. In addition, unlike 
laser glasses in which femtosecond laser writing is easily used to 



produce 3D waveguiding structures, in crystalline laser materials it 
is considerably difficult since the femtosecond laser irradiation typ- 
ically produces negative index decreases in the tracks. Our work 
suggests a solution to realize lasing with 3D beam manipulation in 
crystalline platforms. As a simple example, our simulation has shown 
that, by using 3-Element integration, a 3D 1X4 beam splitter can be 
fabricated under the similar approach (see Supplementary Fig. 3). 
Therefore, one could expect highly compressed laser devices for light 
guiding and beam manipulation, to satisfy some special require- 
ments for particular applications. On the other side, this technique 
could be extended to other doped YAG crystals for ultra-broadband 
applications. For example, if the platforms of Tm:YAG or Er:YAG 
are used 32,33 , the generated laser wavelength could be extended to 
mid-infrared. In case of Yb:YAG chips 34 , the mode-locking of the 
waveguide laser enables the generation of ultrashort pulses of fem- 
tosecond temporal durations. These new miniature laser sources 
with different wavelength or operation regimes could be utilized in 
broader-scope applications of both scientific researches and human 
life. Finally, one could exploit the application of this technique to 
other crystals to obtain waveguide devices, such as frequency con- 
verters (second harmonic generation), frequency comb generators, 
integrated elements for biophotonics, etc. In conclusion, this work 
demonstrates a novel family of monolithic photonic structures and 
their performance as beam shapers or as laser sources with controlled 



SCIENTIFIC REPORTS | 4:5988 | DOI: 10.1038/srep05988 



5 



output properties. Strong capability of fabrication of various photo- 
nic devices by this simple and flexible technique is also shown in a 
designable manner for circuit production. 

Methods 

Fabrication of photonic structures. To fabricate the photonic structures in the 
Nd:YAG crystal we use the ultrafast- laser direct inscription with the sample- 
scanning approach. The Nd:YAG (doped by 1 at.% Nd 3+ ions) crystal wafer is cut 
with dimension of 9.9 X 10 X 2 mm 3 and optically polished. An amplified 
Ti:Sapphire femtosecond laser (Spitfire, Spectra Physics) that delivered linearly- 
polarized pulses with a temporal duration of 120 fs, a central wavelength of 
795 nm and operating at a repetition rate of 1 kHz, is used as the laser source. The 
beam is focused by a 20 X microscope objective and the average power is reduced 
to 1.2 mW (measured after the objective) by using a set of half-wave plate and 
linear polarizing cube, and a calibrated neutral density filter. The laser irradiation 
is controlled with a mechanical shutter. The sample is placed on a XYZ micro- 
positioning stage that allowed scanning the sample at constant velocity (500 um/s) 
while irradiating from the surface of 9.9 X 10 mm 2 with the femtosecond pulses at 
certain depth beneath the surface (—150 um): under our experimental conditions, 
a damage track is then produced with a transverse length of —10 um. A number 
of parallel scans of the sample are done at different heights in order to obtain the 
desired hexagonal photonic structures. 

Characterization of guidance and fluorescence. The waveguiding performance is 
carried out under a typical end-face arrangement. A linearly polarized CW solid state 
laser with single wavelength of 1064-nm is utilized as the light source. The polarized 
light beam is focused by a 20 X microscope object lens, and the coupling into the 
photonic structures. Another microscope object lens is used as the out-coupler, 
through which the transmitted light is collected and imaged by a CCD camera. The 
light powers are measured by a few powermeters. The propagation losses of the 
photonic structures are determined by directly measurement of the light powers from 
the input and output end-faces. 

The u-PL spectra are excited by a 488-nm CW diode laser. The light beam is 
focused in the sample by a 50X microscope objective (N.A. —0.55), thus yielding a 
spot with a diameter of approximately 0.53 um. The subsequent emission of Nd 3+ 
ions is collected by the same microscope objective and, after passing through filters, 
lenses and pinholes, is analyzed by a high- resolution spectrometer. The software 
LabSpec® and WSxM® are used to analyze and fit the spectra and process the 
fluorescence images in terms of the spectral intensity, energy shift, and width of the 
fluorescence lines 35 . 

Refractive index reconstruction. The refractive index changes of the structures are 
obtained by the procedure described in Ref. 27. Briefly, the refractive index map is 
constituted by the superposition of the damage induced refractive index modification 
(reduction) and the stress field induced refractive index map (refractive index 
increment at compressed volumes). The first is obtained from the fluorescence images 
based on the fluorescence intensity as the fluorescence intensity reduction in respect 
to the bulk material is, in a first order approximation, assumed proportional to the 
laser-induced damage in the Nd:YAG network. The conversion of laser induced- 
damage into refractive index reduction is performed by a reference Nd:YAG sample 
that contains different damage tracks fabricated under the same conditions of 
femtosecond laser irradiation. The stress -induced refractive index map has been 
calculated from the fluorescence map obtained in terms of the induced spectral shift 
of Nd 3+ fluorescence lines. From the spectral shift, one calculate the spatial variation 
of residual stress and from this the relative change in the unit cell volume (i.e. the 
lattice dilatation/densification). Finally, the refractive index modification is calculated 
by applying the Claussius-Mossitti approximation 27 . 

Lasing generation. The laser operation experiments in CW and pulsed regimes are 
both performed by utilizing an end pumping system. A polarized light beam at a 
wavelength of 808 nm is generated from a tunable CW Ti:Sapphire laser 
(Coherent MBR-100). A convex lens with focus length of 25 mm is used to couple 
the laser beam into the photonic structures; the radius of the incident pump laser 
beam is approximately 0.75 mm at 1/e 2 position. For the CW waveguide laser 
systems, an input mirror (with the transmission of 98% at 808 nm and the 
reflectivity >99% at —1064 nm) and an output mirror (with the reflectivity >99% 
at 808 nm and 60% at —1064 nm, respectively), are adhered to the input and 
output end face of the waveguide respectively, constructing the Fabry- Perot lasing 
resonant cavity. Here, a set of mechanical bindings are employed to make the SA 
closely attached to the output facet of the waveguide sample. For the Q-switched 
pulsed laser generation, the output mirror is replaced by a bilayer graphene thin 
film coated on a fused silica glass plate as the SA. The graphene film is grown by 
chemical vapor deposition (CVD) on copper and nickel firstly, and then 
transferred to a fused silica glass. The transmission of the graphene film coated 
output mirror is measured to be —70% at 1064 nm. The generated lasers are 
collected by utilizing a 20X microscope objective lens (numerical aperture N.A. — 
0.4) and imaged by using an infrared CCD. We use a spectrometer with resolution 
of 0.2 nm to analyze the emission spectra of the laser beams. The schematic plot 
of the experimental setup for the CW and pulsed laser oscillations are depicted in 
Fig. 4a. 
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